Mon. Not. R. Astron. Soc. OOP. [Tlf23l f 20081 Printed 4 May 2009 (MN I^TeX style file v2.2) 



On the Origin of Exponential Galaxy Disks 



Aaron A. Button 

^ UCO/Lick Observatory, University of California, Santa Cruz, CA 95064 



submitted to MNRAS 



ABSTRACT 

One of the most important unresolved issues for galaxy formation theory is to under- 
stand the origin of exponential galaxy disks. We use a disk galaxy evolution model 
to investigate whether galaxies with exponential surface brightness profiles can be 
produced in a cosmologically motivated framework for disk galaxy formation. Our 
model follows the accretion, cooling, and ejection of baryonic mass, as a function of 
radius, inside growing dark matter haloes. The surface density profile of the disk is 
determined by detailed angular momentum conservation, starting from the distribu- 
tion of specific angular momentum as found in cosmological simulations. Exponential 
and quasi-exponential disks can be produced by our model through a combination of 
supernova driven galactic outflows (which preferentially remove low angular momen- 
tum material), intrinsic variation in the angular momentum distribution of the halo 
gas, and the inefficiency of star formation at large radii. We use observations from 
the SDSS NYU-VAGC to show that the median Sersic index of late-type galaxies is 
a strong function of stellar mass. For blue galaxies, low mass galaxies have n ~ 1.3, 
while high mass galaxies have n ~ 4, with a transition mass of Mgtar — 2.5 x 10^^ Mq. 
Our model with energy driven outflows correctly reproduces this trend, whereas our 
models with momentum driven outflows and no outflows over predict the Sersic in- 
dices in low mass galaxies. We show that the observed fraction of "bulge-less" ex- 
ponential galaxies is a strong function of stellar mass. For Milky- Way mass galaxies 
{Vrot — 220kms~^,Mstar — IO^^Mq) less than 0.1% of blue galaxies are bulge-less, 
whereas for M33 mass galaxies (K-ot — 120 kms~^, Mgtar — lO^^M© bulge-less and 
quasi-bulge-less galaxies are more common, with 45% of blue galaxies having Sersic 
index n < 1.5. These results suggest that the difficulty of hierarchical formation mod- 
els to produce bulge-less Milky- Way mass galaxies is, in fact, not a problem. However, 
the problem of producing M33 like galaxies remains, and will provide a key test for 
hierarchical galaxy formation models, and feedback models inparticular. 

Key words: galaxies: formation - galaxies: fundamental parameters - galaxies: haloes 
- galaxies: kinematics and dynamics - galaxies: spiral - galaxies: structure 



In the current paradigm for galaxy formation, set forth 
by White & Rees (1978) and Fall & Efstathiou (1980), 
disk galaxies are considered to form through the cooling 
of baryons inside dark matter haloes that grow by means 
of gravitational instability and acquire angular momentum 
from cosmological torques. In this paradigm the baryons 
acquire the same distribution of specific angular momen- 
tum as the dark matter, and this is conserved during cool- 
ing. Analytic models based on these assumptions have been 
able to produce disk sizes that are in reasonable agreement 
with observations (e.g. Dalcanton, Spergel & Summers 1997; 



1 INTRODUCTION 



Mo, Mao & White 1998; de Jong & Lacey 2000; Button 
et al. 2007). In this picture the large scatter in disk sizes 
(including the existence of low surface brightness galaxies) 
is a natural consequence of the large scatter in halo spin 
parameter. 



* duttonOucolick.org 



However, detailed hydro-dynamical simulations aimed 
at investigating the process of galaxy formation have in- 
dicated a potentially important problem for the cold dark 
matter (CDM) model. In CDM cosmologies, haloes form hi- 
erarchically by the merging of many lower mass haloes. Be- 
cause the cooling in dense, low-mass haloes is very efficient, 
the baryons in these systems have already cooled by the time 
they merge with the more massive protogalaxy. They reach 
the center of the potential well by means of dynamical fric- 
tion, through which they lose a significant fraction of their 
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angular momentum to the dark matter. Consequently, the 
disks that form in these simulations are an order of mag- 
nitude too small (Navarro & Benz 1991; Navarro & White 
1994; Navarro & Steinmetz 2000). This is known as angular 
momentum catastrophe. 

Attempted solutions to this problem include stellar 
feedback and/or ionizing background radiation which can 
prevent the cooling of the gas in low-mass haloes, thereby 
reducing the angular momentum loss (e.g. Domfnguez- 
Tenreiro, Tissera & Saiz 1998; Weil, Eke & Efstathiou 1998; 
Sommer-Larsen, Gelato & Vedel 1999; Eke, Efstathiou & 
Wright 2000; Mailer & Dekel 2002). An ahernative sug- 
gestion has been to reduce the number of low mass haloes 
by reducing the amplitude of the power spectrum on small 
scales (e.g. Kamionkowski & Liddle 2000; Sommer-Larsen 
& Dolgov 2001). However, these effects have only partially 
successful in resolving the angular momentum catastrophe. 
More recent numerical simulations have shown that in ad- 
dition to feedback, very high numerical resolution is also re- 
quired in order to avoid numerical angular momentum losses 
(e.g. Kaufmann et al. 2007; Ceverino & Klypin 2007; Mayer 
et al. 2008). 

A second problem, related to angular momentum, is to 
understand the exponential density profile of stellar galac- 
tic disks. If disk galaxies form as envisaged in the stan- 
dard picture their resulting surface density distribution is 
directly related to the specific angular momentum distribu- 
tion of the protogalaxy. Motivated by the work of Mestel 
(1963), Crampin & Hoyle (1964) and Gunn (1982), Dalcan- 
ton et al. (1997) made the assumption that the protogalaxy 
has the AMD of a uniform sphere in solid body rotation, 
and showed that the resulting disks are more concentrated 
than exponential. 

Firmani & Avila-Reese (2000) and van den Bosch 
(2001) refined this model by following the build-up and evo- 
lution of the disk density profile within a (cosmologically 
motivated) growing dark matter halo. In these models each 
shell of collapsing material is assumed to have the angular 
momentum distribution of a shell in solid body rotation, and 
the normalization of the angular momentum distribution is 
calculated by assuming that the spin parameter of the halo 
is independent of redshift. In the Firmani & Avila-Reese 
(2000) model the baryonic mass fraction of the galaxy was 
treated as a free parameter, with a fiducial value of 0.05. 
In the standard CDM cosmology (i.e. flm.o ~ 1), such a 
baryon fraction is consistent with the universal value. How- 
ever, with the currently favored a ACDM cosmology (where 
f2ni,o — 0.26) the baryon fraction is ~ 0.17, and the low 
baryons fractions in observed galaxies require an explana- 
tion. In the van den Bosch (2001) model the galaxy mass 
fraction is determined by the efficiencies of cooling and su- 
pernova (SN) driven outflows. Cooling is very efficient in 
galaxy mass halos, thus in order to produce the low galaxy 
mass fractions that are observed (e.g. Hoekstra et al. 2005; 
Mandelbaum et al. 2006) and required to reconcile the stel- 
lar and halo mass functions (e.g. Yang et al. 2007; Conroy 
& Wechsler 2008) galactic outflows are needed. 

The van den Bosch (2001) model highlights the prob- 
lem of producing pure exponential stellar galactic disks in 
a cosmological framework. In these models the disk may 
be approximately exponential over a few scale lengths, but 
the disk is invariably more concentrated than exponential 



in the center. A potential short-coming of the Firmani & 
Avila-Reese (2000) and van den Bosch (2001) models is the 
assumption that the accreting mass shells have the AMDs 
of shells in solid body rotation. 

Bullock et al. (2001b) determined the distribution of 
specific angular momentum in ACDM haloes from cosmo- 
logical N-body simulations, and again concluded that these 
distributions will form overly concentrated disks. This excess 
of low angular momentum material has been confirmed by 
a number of subsequent studies (van den Bosch et al. 2002; 
Chen et al. 2003; Sharma & Steinmetz 2005). These simu- 
lations also indicated that CDM haloes have an excess of 
high specific angular momentum material, relative to that 
required to make an exponential disk. Furthermore, not even 
a single bulge-less disk galaxy with an exponential stellar 
density profile has been formed in hydrodynamical simula- 
tions in the ACDM paradigm. Thus one of the most impor- 
tant unresolved issues in galaxy formation is to understand 
the origin of the exponential proflle of stellar galactic disks. 

There are several reasons why the speciflc angular mo- 
mentum profile of the stellar disk may be different than that 
of the dark matter and hot gas. Firstly only a fraction of the 
hot gas will be able to cool; secondly some of the cooled gas 
may be ejected from the galaxy or progenitor galaxies (e.g. 
Mailer & Dekel 2002) via SN driven winds; thirdly only a 
fraction of the cooled gas will be converted into stars. Fur- 
thermore the surface brightness profile in optical light may 
not reliably trace the surface density profile in stellar mass. 
Thus an important question is whether the processes of star 
formation, feedback and stellar populations can produce ex- 
ponential surface brightness profiles. If this is not possible 
then it signals that the cold baryons acquire a different dis- 
tribution of specific angular momentum than the dark mat- 
ter, or more fundamentally signals a failure of the CDM 
paradigm of structure formation. 

In this paper we use the models described in Dutton & 
van den Bosch (2008), which are an updated version of the 
models used by van den Bosch (2001, 2002) to investigate 
the origin of the exponential light profiles of galactic disks. 
This model includes the detailed cooling, star formation and 
ejection of g function of radius, inside growing dark 

matter haloes with cosmologically motivated density profiles 
and angular momentum distributions. 

This paper is organized as follows. In §2 we give an 
overview of the models, in §3 we discuss the impact of the 
initial angular momentum distribution, star formation, feed- 
back, halo mass, adiabatic contraction and stellar popula- 
tions on the structure of galaxy disks. In §4 we use observa- 
tions from the SDSS to determine how frequent bulge-less 
galaxy disks are, and in §5 we make a comparison between 
these observations and Monte Carlo model samples. We dis- 
cuss our results in §6 and give a summary in §7. 



2 DISK GALAXY EVOLUTION MODELS 

Here we give a brief overview of the disk galaxy evolution 
model used in this paper. This model is described in detail 
in Dutton & van den Bosch (2008; hereafter DB08). The 
key difference with almost all disk evolution models is that 
in this model the inflow (due to gas cooling), outflow (due to 
SN driven winds), star formation rates, and metallicity, are 
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computed as a function of radius, rather than being treated 
as global parameters. The main assumptions that character- 
ize the framework of our models are the following: 

(i) Dark matter haloes around disk galaxies grow by the 
smooth accretion of mass which we model with the Wechsler 
et al. (2002) mass accretion history (MAH). The shape of 
this MAH is specified by the concentration of the halo at 
redshift zero; 

(ii) The structure of the halo is given by the NFW pro- 
file (Navarro, Frenk, & White 1997), which is specified by 
two parameters: the mass and concentration. The evolu- 
tion of the concentration parameter is given by the Bullock 
et al. (2001a) model with parameters (Maccio et al. 2008) 
for a WMAP 5th year cosmology (Dunkley et al. 2009); 

(iii) Gas that enters the halo is shock heated to the virial 
temperature, and acquires the same distribution of specific 
angular momentum as the dark matter. We use the angu- 
lar momentum distributions of the halo as parametrized by 
Bullock et al. (2001b) and Sharma & Steinmetz (2005); 

(iv) Gas cools radiatively, conserving its specific angular 
momentum, and forms a disk in centrifugal equilibrium; 

(v) Star formation occurs according to a Schmidt (1959) 
type law on the dense molecular gas, which is computed 
following Blitz & Rosolowsky (2006); 

(vi) SN feedback re-heats some of the cold gas, ejecting 
it from the disk and halo; 

(vii) Stars eject metals into the inter stellar medium, en- 
riching the cold gas. 

(viii) Bruzual & Chariot (2003) stellar population synthe- 
sis models are convolved with the star formation histories 
and metallicities to derive surface brightness. 

Each model galaxy consists of five mass components: dark 
matter, hot halo gas, disk mass in stars, disk mass in cold 
gas, and ejected gas. The dark matter and the hot gas are 
assumed to be distributed in spherical shells, the disk mass 
is assumed to be in infinitesimally thin annuli. Throughout 
this paper we refer to R as radius, t as time (where t = is 
defined as the Big Bang) and z as redshift. 

For each galaxy we set up a radial grid with 200 bins 
quadratically sampled from between 0.001 and 1 times the 
redshift zero virial radius. As described in DB08 for the pur- 
pose of conserving angular momentum it is more convenient 
to use a grid in specific angular momentum, j. We convert 
the initial grid in R, to a grid in j using j^/G — RA4{R), 
where M{R) is the halo mass enclosed within a sphere of ra- 
dius R. We follow the formation and evolution of the galaxy 
using 400 redshift steps, quadratically sampled from z — 10 
to z = 0. For each time step we compute the changes in the 
various mass components in each radial bin. 

2.1 Angular Momentum Distribution 

The radius at which the cooled gas ends up depends on its 
specific angular momentum, j. Van den Bosch (2001; 2002) 
assumed the j-distribution to be that of a shell in solid body 
rotation. In this paper we adopt cosmologically motivated 
specific angular momentum distributions (AMD). 

The AMD can be specified by 2 parameters: a normal- 
ization (A) and a shape (p, or a). Both the normalization 
and shape parameters are log-normally distributed, with sig- 
nificant scatter. We assume that the spin and shape param- 



eters are uncorrelated, although Bullock et al. (2001b) show 
that there may be a weak correlation between fi and A, which 
needs to be investigated further. 

The spin parameter is defined by 



A = 



t-^vir I -E'vii 



1/2 



GAC; 



5/2 



(1) 



where Mvir, Jvir, and i?vir are the mass, total angular mo- 
mentum and energy of the halo, respectively. 

Bullock et al. (2001b; hereafter BOl) showed that dark 
matter haloes in a ACDM cosmology have specific angular 
momentum profiles that are well fitted by 



M{< j) 

Mvir 



j/jn 



j/jn 



(2) 



Here M(< j) is the halo mass with specific angular mo- 
mentum less than j, Mvir is the halo's virial mass, ^ is a, 
free parameter, and jmax is the maximum specific angular 
momentum in the halo. The value of jmax is given by 



Jmax — 



Jvir/Mv 



(M-l)[-Mln(l-M-^)-l]' 



(3) 



where Jvir/Afvir is determined by the spin parameter. BOl 
found ^ — 1 to be log-normally distributed with a mean 
and standard deviation in log(/i — 1) of ~ —0.6 and ~ 0.4, 
respectively. 

Sharma & Steinmetz (2005, hereafter SS05) used a se- 
ries of non radiative N-body/SPH simulations in a ACDM 
cosmology to study the growth of angular momentum in 
galaxy systems. SS05 showed that Eq. [2] is unable to de- 
scribe the specific angular momentum distribution required 
to make an exponential disk in an NFW halo. They intro- 
duced an alternative function that is able to describe the spe- 
cific angular momentum distribution of exponential disks, as 
well as those of the gas and dark matter in their simulations. 

m(j)=7(a,^), M = 4P^- (4) 

where 7 is the incomplete gamma function. SS05 found the 
distributions of a is log-normally distributed with mean 
log a ~ —0.05, standard deviation in log a ~ 0.11. 



2.2 Conservation of Angular Momentum and 
Halo Contraction 

In order to ensure that specific angular momentum is con- 
served in our model, rather than keeping track of how much 
mass is at a given radius, R, we keep track of how much mass 
is at a given specific angular momentum, j. Thus, under 
the simplifying assumption that Kirc = [GM{< R)/R]^^^, 
where M(< ii) is the total mass within a spherical radius, 
R, the radius that corresponds to a given j is given by 



R = 



J 



GM{< j) 



(5) 



This has a number of desirable properties: 1) At each time 
step it is trivial to calculate how much cold gas is added to 
each bin in j. 2) Over time, as the potential well changes, 
the specific angular momentum of the baryons is automat- 
ically conserved; 3) The adiabatic contraction (Blumenthal 
et al. 1986) of the halo in response to the cooling of the 
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baryons is automatically taken into account. 4) The result- 
ing radial grid is adaptive, as the mapping between j and R 
depends on the amount of mass enclosed. 

As discussed in Button et al. (2007) and DB08, models 
with halo contraction (and standard stellar IMF's) are un- 
able to reproduce the zero points of the V MR relations as 
well Eis the low galaxy formation efficiency required to recon- 
cile the halo mass function and galaxy stellar mass function. 
While there are processes such as dynamical friction and 
impulsive mass loss that can expand the halo, implementing 
these in a galaxy evolution model is a non-trivial task. Thus 
for simplicity we wish to consider a model in which the dark 
halo does not respond to galaxy formation. To calculate the 
mapping between j and radius, for the case of no adiabatic 
contraction, we solve the equation 



R 



Mhalo(-R) +MdiBk(< j) 



(6) 



where Mhaio(J^) is the mass (within a spherical radius R) of 
the dark matter plus hot gas halo in the absence of galaxy 
formation, and Mdisk(< j) is the mass of the disk (gas plus 
stars) with specific angular momentum less than j. This 
way the self-gravity of the disk is included but adiabatic 
contraction is ignored. 



2.3 Star Formation 

Following Blitz & Rosolowski (2006) we assume that star 
formation takes place in dense molecular gas, traced by hy- 
drogen cyanide (IfCN), with a constant star formation effi- 
ciency: 

SsFR eSF Smol.HCN 



[MQpc-2Gyr-i] [Gyr-i] [Mopc" 



(7) 



where ^sf ~ 10 - 13Gyr"^ (Gao & Solomon 2004, Wu 
et al. 2005). Expressing this equation in terms of the total 
gas content: 

SsFR eSF Sgi 



,fmol 72.HCN , 



(8) 



[M0pc-2Gyr-i] [Gyr-i] [Mopc" 

where 7?.hcn = Smoi.HCN/Smoi is the ratio between the 
dense molecular gas (as traced by HCN) and the total molec- 
ular gas, and /moi is the molecular gas fraction. 

The fraction of gas that is molecular, /moi, can be de- 
fined in terms of the mass ratio between molecular and 
atomic gas, 7?.moi by 



/moi — 



(9) 



Empirically Wong & Blitz (2002) and Blitz & Rosolowsky 
(2004; 2006) have argued that T^moi is determined to first 
order by the mid plane pressure, Pcxt. The most recent re- 
lation from Blitz & Rosolowsky (2006) is 



4.3 ± 0.6 X 



(10) 



where k is Boltzmann's constant, and Pcxt/fc is in cgs units. 
For a gas plus stellar disk the mid plane pressure is given, 
to within 10% by (Elmegreen 1993) 



Csta 



(11) 



where cTgas and (Tatar are the velocity dispersions of the gas 
and stellar disk, respectively. For simplicity we will assume 

Tgas /fstar =0.1. 

Based on the arguments and references in Blitz & 
Rosolowski (2006) we adopt the following relation for 7?.hcn 



TlucN = 0.1 1 



[200 Mepc- 



(12) 



In the low pressure, low molecular density regime, 7?.hcn 
0.1, and thus Eq. ([7|) asymptotes to 



jSFR 



esF 



0.013 



[M0pc-2Gyr-i] [Gyr-i] [Mq pc- 



[A/opc- 



.(13) 



In the high pressure, high molecular density regime, 7?.hcn oc 
E^'oi, and eq.((7)l asymptotes to the familiar Schmidt- 
Kennicutt (SK) relation 



jSFR 



esF 



0.012 



[Af0pc-2Gyr-i] [Gyr-i] [Mq pc" 



[Mopc 



.(14) 



Furthermore, with esF = 13 Gyr^^, we recover the coefficient 
of esF = 0.16 of the standard SK relation. 

We implement the star formation recipe given by Eq.([7| 
as follows. At each time step and annulus in the disk, we 
calculate the star formation rate. Then we use the following 
approximation (valid for times steps small compared to the 
star formation time scale) to calculate the mass of newly 
formed stars 

AAfetar(i?) = A{R) Esfr(-R, t) At, (15) 
with A the area of the annulus, and At the time step interval. 



2.4 Supernova Feedback 

We assume that SNe drive mass outflows, and that the out- 
flows move at the local escape velocity of the disk-halo sys- 
tem. This choice is motivated by the fact that it maximizes 
the mass loss from the the disk-halo system (lower velocity 
winds will not escape the halo, and higher velocity winds 
will carry less mass). 

For our energy driven wind model following Dekel & 
Silk (1986) we assume that the kinetic energy of the wind is 
equal to a fraction, eefb , of the kinetic energy produced by 
SN. However, contrary to Dekel & Silk (1986) we apply this 
energy condition locally in the disk as a function of radius, 
rather than globally to the whole galaxy. Thus the mass 
ejected from radius, R, during a given time step is given by 

2 eefb Esn J?sn 



AMejcct(i?) 



Klc(i?) 



-AM,t.r{R). 



(16) 



Here AMstar(-R) is the mass in stars formed at radius, R, 



Esn 



10^^ erg 



5.0 X lO'^ km^ s~^ Mq is the energy pro- 
duced by one SN, and ?7sn = 8.3 x 10"'^ is the number of SN 
per solar mass of stars formed (for a Chabrier (2003) IMF). 

Similarly for our momentum driven wind model we as- 
sume that the momentum of the wind is equal to a fraction, 
fMFB, of the momentum produced by SN, thus the mass 
ejected from radius, R, during a given time step is given by 



AMejcct(i?) 



fMFB PSN r]SN 



AMetar(P). 



(17) 



Kec(i?) 

Here psN = 3 x W^Mq kms~^ is the momentum produced 
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by one SN, assuming that each SN produces ~ IOMq of 
material moving at t; ~ 3000 kms~^ (Murray, Quataert & 
Thompson 2005). Note that this corresponds to a kinetic 
energy of 4.5 x lO'' M© km^ s"^ ~ 10^^ erg. 

We assume that the ejected mass is lost forever from the 
system: the ejected mass is not considered for later infall, 
and the corresponding metals are not used to enrich the 
infalling gas. This is clearly a dramatic oversimplification, 
but we make this choice to maximize the amount of gas 
that is lost from the system. 

2.5 Overview of Parameters 

The input parameters of our models are as follows. 

(1) Cosmology: fini.cf^A, S7b, o"8, h, n. In this paper 
we adopt a flat ACDM cosmology motivated by the 5th 
year WMAP results (Dunkley et al. 2009), with fJ^.o = 
0.258, t^A = 0.742, f7b = 0.044, ag = 0.80, ft = 0.7, and 
n = l. 

(2) Halo structure: K,F,ainc- We adopt the Bullock 
et al. (2001a) model with F = 0.01, K = 3.7, and cri„c = 
0.25. These parameters reproduce the distribution of halo 
concentrations of relaxed dark matter haloes in cosmolog- 
ical N-body simulations (e.g. Wechsler et al. 2002; Maccio 
et al. 2007; 2008) 

(3) Angular momentum distribution: A, (Jin a, ct (or fi), 
(or criog(,j-i)). As fiducial values we adopt a median 

spin parameter A = 0.025 with a scatter crinA ~ 0.35 (see 
DB08). For the angular momentum distribution we consider 
two fitting functions: the Bullock et al. (2001b) formula 
with median fi — 1 — 0.3, criog(n-i) = 0.45; and the Sharma 
& Steinmetz (2005) formula with a median a = 0.9 and 
o"ioga = 0.11. For our fiducial model we adopt the Sharma 
& Steinmetz (2005) angular momentum profile. 

(4) Star formation: esp. We use a star formation model 
based on dense molecular gas, which we compute following 
Blitz & Rosolowsky (2006). We adopt ?sf = 13Gyr-\ 

(5) Feedback: eEFB, cmfb, Vsn, -Esn, psn- We adopt an 
energy per SN of Esn = 10^^ erg, a momentum per SN of 
PSN = 3 X 10* M© kms~\ a SN rate of 77SN = 4 x 10"^ per 
solar mass of stars formed. We treat eBPB and emfb as free 
parameters, with fiducial values of eefb ~ 0.25, emfb ~ 0. 

(6) Stellar populations and chemical enrichment: 
TZ,y, Z\iat, and the choice of initial mass function (IMF). 
We adopt the Chabrier IMF, a return fraction TZ = 0.35, 
and a stellar yield y = 0.02, and a hot gas metallicity of 
0.002 (i.e. Q.IZq). 



3 SURFACE DENSITY PROFILES 

The surface density profile of a centrifugally supported disk 
is determined by two things: 

(i) the distribution of specific angular momentum 
(AMD), j, of the baryons (i.e. stars and cold gas); and 

(ii) the shape and amplitude of the rotation curve. 

The AMD of the baryons may be modified from that pro- 
duced by tidal torques, by inefficient cooling and SN driven 
outflows. The difi'erence between the surface density pro- 
file of the stellar disk and baryonic disk is then determined 
by the (radially dependent) efficiency at which the cold 



gas is turned into stars. Another complication is that ob- 
served light is a biased tracer of the stellar (and baryonic) 
mass. Below we discuss each of these effects for a fidu- 
cial model galaxy. This model has a total virial mass of 
Mvir = 9.0 X 10^^ M0, a median halo concentration of this 
mass, and unless other- wise specified a halo spin parame- 
ter of A = 0.025. This model also has no halo contraction, 
as models with halo contraction are unable to reproduce 
the zero points of the TuUy-Fisher (TF) (TuUy & Fisher 
1977) and size- mass relations (DB08). However, as we show 
in § 333] &[5]4] none of our main results are sensitive to this 
assumption. 

3.1 The Role of the Initial Angular Momentum 
Distribution 

In the standard picture of disk galaxy formation the dis- 
tribution of j of the cooled gas is the same as that of the 
hot gas and dark matter. As discussed in § [51 N-body sim- 
ulations have shown that there is a universal j profile, de- 
scribed by one parameter in BOl; a in SS05). In Fig[l] 
we show the effect of the initial AMD shape parameter on 
the surface density and circular velocity profiles for our fidu- 
cial model galaxy. The upper panels use the BOl AMD with 
^ — 1.04, 1.11, 1.30, 1.90, and 3.50, while the lower panels 
use the SS05 AMD with a = 0.54, 0.69, 0.90, 1.15, and 1.48. 
For both AMD's these values correspond to the —2, —1, 0, 1, 
and 2(7 of the distribution found in cosmological simulations. 

In order to isolate the effects of the AMD from the 
effects of feedback, we fix the galaxy mass fraction and 
galaxy spin parameter to be equal to the values for a galaxy 
with a — 0.9 and our favored energy driven feedback model 
(eEFB = 0.25): nig^i = 0.058 and Agai = A = 0.028. 

Both the BOl and SS05 AMD's resuh in roughly the 
same surface density profiles and rotation curve shapes. The 
generic result is for lower values of fi and a to result in 
higher surface densities and circular velocities at small radii. 
As /i and a increase the central densities decrease and the 
circular velocity profiles become flatter. There are however 
two notable differences. The BOl profile has a maximum j, 
which results in a truncation in the surface density profile, 
this truncation occurs at a smaller radius for larger values of 
jj,. The SS05 profile has less low angular momentum material 
for high values of a than the BOl angular momentum profile 
for high values of /i. 

SS05 showed that the AMD at large j depends on the 
method used to calculate the AMD. In particular, the cell 
method used by BOl tends to result in a truncation of the 
AMD at large j, which is better fitted by the BOl profile. 
The particle method used by SS05 results in a smooth and 
extended AMD at large j, which is better fitted by the SS05 
profile. If one removes the outer 5% of the j'-distribution, 
SS05 showed that both profiles provide equally good fits to 
the AMDs computed using the cell and particle methods. 

The concentration parameters for each model are given 
in Fig. [1] The concentration parameter of the stellar (C82s) 
and baryonic disks (C82t) are defined as the ratio between 
the radii enclosing 80 and 20% of the disk stellar/baryonic 
mass. For reference an exponential profile (Sersic n — 1) 
corresponds to a concentration of 6*82 = 3.6, a Sersic n = 2 
profile corresponds to a C82 ~ 5.8, and a de Vaucouluers 
profile (Sersic n = 4) corresponds to a Cs2 ~ 11.5. 
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Figure 1. Effect of angular momentum shape parameter on surface density and circular velocity profiles. Radii are expressed in terms 
of the virial radius, R^ii, and velocities in terms of the circular velocity at the virial radius, V^ir- For the cosmology we adopt, the virial 
radius (in units of h^^kpc) is ~ 2.068 times the virial velocity (in units of kms~^). All models have mg^i = 0.058, A = 0.028, and 
Mvir = 10"-*/i-l Mq. The upper panels use the BOl angular momentum profile (specified by /i), while the lower panels use the SS05 
angular momentum profile (specified by a). The five values of the angular momentum shape parameter correspond to the mean, itlcr 
and ±2(7 of the distributions found in cosmological simulations. The stellar disk is given by the blue solid lines, the total disk by the blue 
dotted lines. The gas disk is given the the solid green lines, the molecular gas is given by the dotted green lines, and the atomic gas is 
given by the dashed green lines. The circular velocity of the halo is given by the solid red line and the total circular velocity is given by 
the solid black line. The circles show the circular velocity at three characteristic radii: 2.15 stellar disk scale lengths (black open circles); 
the radius enclosing 80% of the stellar mass, iJgOs (blue filled circles); and the radius enclosing 80% of the cold gas mass, -RsOc (green 
filled circles). For each model the cold gas fraction, /gas, and stellar disk concentration, C82s, and the baryonic disk concentration, C's2t, 
are given. 



The BOl AMD results in less dispersion in the con- 
centrations of the stellar and baryonic disks than the SS05 
AMD, both at high and low values of a or jj.. It should be 
noted that some of the very high concentrations Cs2s — 100 
produced with the SS05 AMD are unrealistic. These are 
caused by the models having too much mass with low specific 
angular momentum, and our requirement that this specific 
angular momentum be conserved. In order for mass with 
low specific angular momentum to be in centrifugal equilib- 
rium it has to be at small radii. However, when the baryons 
dominate the potential at small radii, this increases the cir- 
cular velocity, which in turn requires the baryons move to 
even small radii, and so on. When centrifugal equilibrium 



is reached, the inner disk is very compact. In practice the 
disk will become unstable before it reaches centrifugal equi- 
librium, and the baryons in the inner part of the galaxy 
will become at least partially supported by random motions, 
which will result in lower concentrations. 



3.2 The Role of Star Formation 

Fig. [1] also shows that the surface density profile of the stel- 
lar disk is very different than that of the total baryonic disk. 
The stellar disks are closer to exponential and have smaller 
sizes than the baryonic disks. These differences are due to 
the increased efficiency of star formation at higher gas den- 
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Figure 2. Effect of feedback on tfie surface density and rotation curves using the SS05 angular momentum profile. The upper panels 
have a = 0.9 with cefb varying from to 1, while the lower panels have eEFB = 0.25 with a varying over the 2a range. The line 
types and colors are the same as in Fig[T] Without feedback (upper panels, far left) the model has too many baryons, which results in 
a disk that is too small, too concentrated, gas poor, and which rotates too fast. As the efficiency of feedback in increase (upper panels, 
from left to right), the baryon fraction decreases, which results in a larger disk, with lower concentration, higher gas fraction and lower 
rotation velocity. A model with maximal feedback (upper panels, far right) has a large, low concentration, gas rich disk, with a rising 
inner rotation curve. 



sities (the Schmidt-Kennicutt law). More specifically, in the 
star formation model that we adopt, stars only form out of 
molecular gas, and the molecular fraction is a strong func- 
tion of the surface density of stars and gas. In Fig. [1] the 
surface density of the molecular gas is given by the dotted 
green lines, and that for the atomic gas by green dashed 
Unes. At small radii the gas is mostly molecular, whereas 
at large radii the gas is mostly atomic, and hence results in 
very inefficient star formation at low gas densities. 

Our models thus predict that gas disks should have 
larger sizes than stellar disks, and that the baryonic disks 
are in general not exponential. Thus a study of the baryonic 
surface density profiles of galaxy disks would provide useful 
constraints to the initial AMDs and the efficiency of star 
formation at low gas densities. 

An interesting consequence of the inefficiency of star 
formation at low gas densities, is that this results in a outer 
disk with a shorter scale length than the main disk. Such 
"disk breaks" are often seen in observed spiral galaxies (e.g. 



Pohlen & Trujillo 2006). Why some galaxies show breaks, 
and others are purely exponential to large radii, and still 
others show anti-truncations (i.e. larger outer disk scale 
lengths), is an open question. If the efficiency of star forma- 
tion at low gas densities varies between galaxies this could 
provide an explanation for the wide range in outer disk pro- 
files observed in spiral galaxies. A study of baryonic disk 
density profiles would thus provide a means of distinguish- 
ing between a star formation efficiency origin vs a dynami- 
cal origin (e.g. Debattista et al. 2006; Younger et al. 2007; 
Roskar et al. 2008; Foyle, Courteau & Thacker 2008) for 
these features. 



3.3 The Role of Feedback 

The models in Fig. [1] had the galaxy mass fraction and spin 
parameter fixed. However, in a halo with Mvir — 9x 10^^ Mq, 
the cooling efficiency is high, such that in the absence of 
feedback the galaxy mass fraction is ~ 80%. In order to 
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Figure 3. Effect of halo mass on surface density and circular velocity profiles in a model with feedback. All models have eefb = 0.25, 
ct = 0.90 and A = 0.025. The redshift zero virial mass (in h = \ units) is given above each column, the corresponding virial velocity 
is given in the velocity panels. The upper panels are for a model without adiabatic contraction, the lower panels include adiabatic 
contraction. The line types and colors are as in Fig. [l] Low mass haloes result in low galaxy mass fractions (mgai), high galaxy spin 
parameters (Ag^i), high gas fractions (/gas), low stellar concentrations (C82s), and dark matter dominated inner rotation curves. As halo 
mass increases rrigai increases (up until Mvir ~ lO^^h"^ Mq), Ag^i decreases, /gas decreases and Cg2s increases, and the inner rotation 
curve becomes dominated by baryons. Adiabatic contraction increases the density of the halo, especially in the inner few percent of the 
virial radius. The differences between models with and without adiabatic contraction are in general small, with the exception of the 
amplitude of the inner rotation curve. See text for further details. 



produce galaxies with realistic mass fractions < 33%, mass 
outflows (i.e. feedback) are required. 

The effect of feedback on the surface density profiles is 
shown in Fig[2l The upper panels show a model with the 
SS05 AMD and q = 0.90 and with the feedback efficiency 
varying from left to right: eefb = 0.00,0.10,0.25,0.50, and 
1.00. The lower panels show a model with eefb = 0.25, and 
the SS05 AMD varied over the 2a range as in Fig. [l] (this 
corresponds to model III in DB08). 

Feedback effects the surface density profile in two ways: 

1) it removes mass, decreasing the galaxy mass fraction, and 

2) it preferentially removes low angular momentum mate- 
rial which also modifies the AMD. The galaxy mass frac- 
tions, mgai = Mgai/Mvir (defined as the ratio between the 
galaxy mass and the total mass of the galaxy-halo system), 
and galaxy spin parameters, Agai = (jgai/?ngai)A, (where 



Jgai = J^ax/Jvii IS the ratio between the total angular mo- 
mentum of the galaxy and that of the galaxy-halo system) 
are given for each model. As the feedback efficiency is in- 
creased the galaxy mass fraction decreases and the galaxy 
spin parameter increases (upper panels of Fig. [2]). Both of 
these result in a more extended disk and a more gradu- 
ally rising rotation curve. By decreasing the amplitude of 
the rotation curve at small radii gas will settle in centrifu- 
gal equilibrium at a larger radius, hence resulting in a less 
concentrated disk. Feedback also preferentially removes low 
angular momentum material, which also directly results in 
a less concentrated disk. 

In these models the halo was uncontracted (as mod- 
els with halo contraction over-predict the zero point of the 
TF relation). Halo contraction would make the trends with 
feedback efficiency stronger, as high galaxy mass fractions 
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will result in more halo contraction, higher circular velocities 
and hence higher stellar central densities. 

The amount of gas ejected as a function of radius de- 
pends on both the star formation rate and the depth of the 
potential well. Since both star formation is more efficient and 
the potential well is deeper in the centers of galaxies it is not 
obvious which effect should dominate. Indeed, as shown in 
DB08, globally, the depth of the potential well dominates 
over the star formation efficiency, such that galaxies in less 
massive haloes eject a higher fraction of their accreted gas. 
However, locally, for a fixed halo mass, smaller (i.e. lower 
spin parameter) galaxies are more efficient at removing their 
baryons, because the higher star formation rates are more 
than enough to compensate for the deeper potential well. 

The lower panels of Fig. [5] show the effect of varia- 
tion in the AMD shape parameter, a, over the 2a range 
as found in cosmological simulations. These models can be 
compared with those in the lower panels of Fig. [T] where 
the galaxy mass fractions and spin parameters were fixed at 
w.gai = 0.058 and Agai = 0.028 respectively. In Fig. [2] the 
galaxy mass fractions and spin parameters are determined 
by the efficiencies of feedback and cooling. Qualitatively the 
surface density profiles and circular velocity profiles are sim- 
ilar between the two models. However, there are some small 
differences which highlight the impact the distribution of the 
baryons have in determining the amount of mass that can be 
removed. Larger values of a result in lower galaxy mass frac- 
tions (i.e. feedback is more efficient). This is a direct result 
of the shallower potential well created by a less concentrated 
baryonic disk, which makes it easier to eject baryons from 
the center of the galaxy. For the model with a = 1.48, a 
comparison between Fig. [T] and Fig. [2] shows that feedback 
has removed the "bulge" component, creating a hollow core 
in the disk density profile. This results in a stellar disk less 
concentrated than exponential. On the other hand, for low 
values of a, there is more low angular momentum material, 
which results in a deeper potential well, from which it is 
harder to remove gas. 

3.4 The Role of Halo Mass 

Having discussed the effects of the AMD and feedback ef- 
ficiency on the surface density profiles for a single halo 
mass, with median concentration and spin parameters, we 
now turn our attention to halo mass dependence of galaxy 
surface density profiles. Fig. [3] shows the surface density 
and circular velocity profiles for 5 model galaxies with 
log(Mvir/[/i"^M0]) = 10.0,11.0,11.5,12.0, and 13.0. AH of 
these models have the same halo angular momentum pa- 
rameters: a — 0.90, A = 0.025, yet their resulting galaxy 
spin parameters, Agai, and stellar concentrations, C82s, vary 
systematically with halo mass. Below a halo mass of Mvir = 
10^^ M0 almost all of the baryons that accrete onto the halo 
have had time to cool, so the differences between the struc- 
tural properties of the galaxies are a result of the efficiency 
at which feedback can remove mass and angular momen- 
tum from the disk and halo. As shown in DB08, mass loss is 
more efficient from lower mass haloes, and thus feedback has 
a larger impact on the density profiles of low mass galaxies. 
Since star formation is biased to small radii (due to the den- 
sity dependence of the star formation rate), the mass that 
gets removed has lower than average specific angular mo- 



mentum. The ejection of low angular momentum material 
directly lowers the concentration of the stars. The reduced 
galaxy mass fractions and higher galaxy specific angular mo- 
mentum result in a more dark matter dominated circular 
velocity profile. 

In massive haloes, A/vir — W^'^h~^ M© feedback is in- 
efficient at removing mass. Cooling is also inefficient, which 
results in a low galaxy mass fraction, nigai = 0.033. But be- 
cause cooling occurs from the inside out, the resulting spin 
parameter of the galaxy is very small, Agaj = 0.008. Thus the 
resulting galaxy is very compact and concentrated, which 
resembles early type galaxies more than late type galaxies. 
However, we do not consider these models to be realistic, 
as the effects of mergers and secular evolution (which are 
not taken into account in our models) are likely to play a 
significant role in producing bulges in high mass haloes. 

3.5 The Role of Adiabatic Contraction 

The lower panels of Fig. |3]show the same models as the up- 
per panels but including the effect of adiabatic contraction. 
Adiabatic contraction increases the density of the dark mat- 
ter halo (red lines in Fig. |3} in the inner 10% of the virial 
radius. A deeper potential well results in a higher escape 
velocity, and thus mass ejection is less efficient. As shown in 
Fig. |3] and as discussed in Dutton & van den Bosch (2008) , 
models with adiabatic contraction have higher rotation ve- 
locities and smaller sizes than models without adiabatic con- 
traction. As might be expected the effect on the density 
profile of the disk is to increase its concentration, but the 
effects are in general small, especially in low mass haloes. 
Thus adiabatic contraction does not play a significant role 
in determining the density profiles of galaxy disks. 

3.6 The Role of Stellar Populations 

Fig.|4]shows surface brightness profiles in K-, I- and B-bands 
of the models in the lower panels of Fig. [2] (varying AMD 
shape) and the upper panels of Fig. |3] (varying halo mass) . 
In each panel the concentration parameter, Cs2, in K-, I-, 
and B-band light as well as in stellar mass is given. Concen- 
trations are largest in stellar mass and smallest in B-band 
light. The differences depend both on the concentration of 
the stars and the mass of the galaxy, with differences typi- 
cally a factor of ~ 1.4. These differences are caused by color 
gradients: inner disks being redder than outer disks. The 
color gradients are a reflection of the inside-out formation of 
stellar disks, both because our star formation law results in 
more efficient star formation at higher gas densities (which 
occur at smaller radii), and because the gas disk is growing 
with time which reflects the growth of the dark matter halo. 
It is interesting to note that there are non-negligible dif- 
ferences between concentrations in K-band light and stellar 
mass, which we return to in ijij 14. 2l fc 15.31 

Fig. |4] also shows the exponential disk flts as dashed 
lines. The scale lengths in stellar mass, K-, I-, and B-band 
light are given for each model. The scale lengths in B- 
band lights are typically 1.14 times larger than those in 
I-band light, and 1.25 times larger than those in K-band 
light. Observationally, for nearby disk galaxies, MacArthur 
et al. (2003) finds a mean ratio of B-band to R-band scale 
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Figure 4. Effect of stellar populations on concentration index and disk scale lengths for model galaxies in the lower panels of Fig. [2] and 
upper panels of Fig. [3] For each model the surface brightness profiles in B-band (blue), I-band (red) and K-band (black) light are shown. 
The dashed lines show the exponential disk fits. The concentration parameters (defined as the ratio between radii enclosing 80 and 20% 
of the light (or stellar mass) and disk scalelengths are given for each model. 



lengths of Rd.s/ Rd,R = 1.15, with a standard deviation of 
0.15, and B-band to H-band scale lengths of i?d,B/^d,H = 
1.34 with a standard deviation of 0.20. Our models are con- 
sistent with these observations. If we also consider that ex- 
tinction will most likely cause the observed gradients to be 
over-estimated, this will bring the models into even better 
agreement with the observations. 

K-band light is generally considered a good tracer of 
stellar mass, but in our models the scale length of the K- 
band hght is typically 1.2 times larger than the scale length 
in stellar mass. This also means that the scale lengths in B- 
band hght are typically 1.5 times larger than those in stel- 
lar mass. Thus in order to derive stellar mass scale lengths 
(which are a fundamental parameter of disk galaxies) from 
observations, it is necessary to take into account the color 
gradients which will result in stellar mass-to-light ratio gra- 
dients, even in the near IR. 



4 HOW COMMON ARE BULGE-LESS 
EXPONENTIAL DISKS? 

We have shown that it is possible to produce quasi- 
exponential disks with ACDM initial conditions. However, 
these galaxies are not typical (in the sense they do not oc- 
cur for median values of the input parameters) for any halo 
mass between 10^^ and 10^^. However, in our models, low 
mass haloes produce disks with closer to exponential stellar 
profiles. Thus, an important observational question is to de- 
termine what fraction of galaxies, and late-type galaxies in 
particular, have bulge-less exponential disks and low bulge 
fractions, as a function of stellar mass. 

It is generally thought that bulge-less disk galaxies are 
fairly common, and as such, this presents a challenge to hi- 
erarchical galaxy formation models. This is especially the 
case if a significant fraction of galaxy bulges were formed 
via secular processes (e.g. Kormendy & Fisher 2005). Such 



a conclusion is supported by recent analyses that find bulge- 
less disk galaxy fractions of 15% for a sample of edge-on 
disk galaxies (Kautsch et al. 2006), and ~ 20% for sam- 
ple of visually classified quasi-bulge-less galaxies (Barazza 
et al. 2008). However, these analyses do not consider any 
mass dependence to the fraction of bulge-less disk galaxies, 
which may lead one to over predict the frequency of massive 
bulge-less galaxies. 

In an attempt to get at least a preliminary (if not a 
definitive) answer to the question of how common are ex- 
ponential bulge-less disks (as a function of stellar mass) we 
make use of Sersic indices from the publicly available low 
redshift New York University Value Added Catalog (NYU- 
VAGC, Blanton et al. 2005). The low-z NYU-VAGC is 
based on the second data release of the SDSS (Abazajian 
et al. 2004), and consists of 28089 galaxies at distances of 
10-200 Mpc (0.0033 < z < 0.05), which have been deter- 
mined by correcting for peculiar velocities. 

The Sersic profile is given by 

E(i?) = Eoexph(i?/i?,)'^"], (18) 

where Eq is the central surface density, Rs is the scale radius, 
and n is the Sersic index. For n — 4 this profile corresponds 
to a de Vaucouleurs profile (de Vaucouleurs 1959), for n = 1 
to an exponential profile, and for n = 0.5 to a gaussian. 
The NYU-VAGC provides Sersic indices measured in each 
of the five SDSS pass-bands: ugriz. The Sersic fits are per- 
formed with circular apertures, and thus are likely to given 
erroneous results for highly inclined disks. Furthermore, the 
surface brightness profiles of highly inclined disks are often 
strongly affected by extinction, so we wish to remove in- 
clined disks remove from our sample. To do this we use the 
axis ratios measured from exponential fits (taken from the 
main SDSS catalog), and remove galaxies with axis ratios 
smaller than b/a = 0.5 (i.e. corresponding to an inclination 
of ~ 60°). This removes about half the galaxies in the sam- 
ple. De Vaucouleurs fit axis ratios are also available, but 
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Figure 5. Color-stellar mass and Serisc index-stellar mass relations for galaxies in the low redsliift NYU-VAGC. The color coding 
corresponds to the relative number of galaxies in each bin, weighted by l/Vmax, on a logarithmic scale. Galaxies are selected to have 
stellar mass Mstar > IO^'^Mq, redshift 0.005 < z < 0.05, and axis ratio b/a > 0.5. We separate galaxies into red (middle) and blue 
(right) according to their position in the u — r color-r magnitude diagram, using the relation in Baldry et al. (2004). In addition we 
require that red galaxies are not bluer than 0.1 dex from the red sequence, and that blue galaxies are at least 0.05 dex bluer than the 
red sequence. We show the median (solid), 15.9th and 84.1th (dashed), 2.3th and 97.7th (dotted) percentiles of the distributions of color 
and Sersic index in stellar mass bins of width 0.2 dex. The error bars show the Poisson error on the median. The thick solid lines show 
fits to the median relations, see Table [T] for fit parameters. For the red galaxies we show the relations for the blue galaxies as thin solid 
lines, and vice versa for the blue galaxies. 



these often give axis ratios more representative of the bulge 
or bar, rather than the disk, so they are not considered useful 
for our purposes. 

For both Sersic index and axis ratio we adopt the i-band 
fits as this provides a trade off between signal to noise and 
wavelength. We calculate stellar masses using the relation 
between g ~ r color and stellar mass-to-light ratio from Bell 
et al. (2003), with an offset of -0.1 dex corresponding to a 
Chabrier (2003) IMF: 

log Metar/L, = -0.406 + 1.097(5 - r) (19) 

The upper panels in Fig. [S] show the relation between 
color and stellar mass. We separate galaxies into red and 
blue galaxies using u—r colors following Baldry et al. (2004), 
who determine the optimal color separation between red and 
blue galaxies to be 

{u ~ r)div = 2.06 - 0.244 tanh (^Ml±^:91^ , (20) 

where Mr is the r-band magnitude of the galaxy. In addition, 
to keep the red-sequence and blue-cloud "pure" in the g — r 
vs Matar plane, we require that red galaxies are not bluer 
than 0.1 dex from the red-sequence, and that blue galaxies 
are at least 0.05 dex bluer than the red-sequence. The color- 



mass relations for red and blue galaxies are shown in the 
upper middle and upper right panels respectively. 

The lower panels in Fig. [5] show the relations between 
Sersic index (measured in the i-band) and stellar mass. For 
all galaxies (left panel) the Sersic index also shows a bi- 
modality, with a cloud of points at high masses around 
n = 4, and at low masses around n — 1.3. There is a 
rapid transition in Sersic index between a stellar mass of 
Msta,! = 1 X 10^° to 5 X 10^" M0. It is tempting to asso- 
ciate this bi-modality with the color bi-modality. However, 
if we look at early and late-types separately, we see the same 
trends: low mass galaxies have low Sersic index and high 
mass galaxies have high Sersic index. We return to the cor- 
relation between color and Sersic index in § 14.21 below. 

In order to calculate residuals, and to ease comparisons 
with models we fit the relations in Fig. [S] with the following. 
We fit the color-mass relation of the red galaxies with a 
single power-law 

(ff-r) = (5-r)o + log(^)". (21) 

We fit the color-mass relation of the blue galaxies with a 
double power-law 
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Figure 6. Fraction of galaxies at a given stellar mass with Sersic index less than 1.0 (cyan), 1.5 (blue), 2.0 (green), 2.5 (red), and 4.0 
(magenta). The left panels show results for all galaxies, middle panels for red galaxies, and the right panels for blue galaxies. The black 
points show the fraction of red galaxies at a given stellar mass. The upper panels show the fractions in linear space. The lower panels 
show the fractions in log-space, which highlights the scarcity of bulge-less galaxies (i.e. low n) at high masses. 



Table 1. Best fit parameters for median color-mass 
and Sersic-mass relations from Fig. \5\ using the func- 
tions in Eqs. I I2U23I I. 
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We fit the Sersic index - stellar mass relations with the fol- 
lowing function: 

log n = log n2 + ^^^o,i„,(M...,/Mo) • (23) 

Here ni is the asymptotic value of n at low mass, n2 is the 
asymptotic value of n at high mass, Mo is the transition 
mass, and 7 controls the sharpness of the transition. The 
best fit parameters are given in table [1] 

From Fig[5]it is apparent that massive galaxies with low 



Sersic indices are rare, both for red and blue sub-samples. 
To be more quantitative Fig. |6] shows the fraction of galax- 
ies at a given stellar mass with Sersic index less than 1.0 
(cyan), 1.5 (blue), 2.0 (green), 2.5 (red), and 4.0 (magenta). 
Also shown is the fraction of all galaxies that are red (black 
dotted). The red-galaxy fraction steadily increases from low 
to high stellar mass. Note that our red galaxy fraction is 
biased high due to our axis ratio cut, which preferentially 
removes disk (and hence typically blue) galaxies. 

At the transition mass of 2.5 X 10^°Mq the median 
Sersic index of the late types is n ~ 2. We adopt a Sersic 
index of n = 1.5 as an upper-limit on the Sersic index of 
a bulge-less galaxy (which was also the definition used by 
Bell 2008). Thus at the transition mass, less than ~ 15% of 
blue galaxies are bulge-less. At higher masses the bulge-less 
galaxy fraction drops rapidly. To show just how rare massive 
bulge- less galaxies are, the lower panels show the fractions 
on a logarithmic scale. 

Using the TF relations from Pizagno et al. (2005) and 
Button et al. (2007), galaxies with rotation velocities of 220 
kms"'^ (e.g. such as the Milky- Way) have a median stellar 
mass of Mstar — IO^^Mq. Note that the actual stellar mass 
of the Milky- Way is thought to be ~ 5 x 10^°Mq (Widrow, 
Pym & Dubinski 2008), which means that the Milky- Way 
does not fall on the TF relation (Hammer et al. 2007). These 
galaxies have a median Sersic index of n ~ 3.5, and less than 
0.1% of these galaxies are bulge-less (i.e. have n < 1.5). 
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Figure 7. Dependence of Sersic index on pass-band. Tlie upper panels sliow the median Sersic index - stellar mass relation for Sersic 
index measured in each of the five SDSS pass-bands: u,g,r,i, and z. The lower panels show the deviation of the Sersic index measured 
in the u,g,r, and z bands relative to that measured in the i band. The dependence of Sersic index on wavelength is strongest for blue 
galaxies, and peaks at the transition mass of Mstar — 2.5 X 10"'^''AfQ. 



This provides some relief for theorists who have been try- 
ing unsuccessfully to produce a bulge-less Milky- Way mass 
galaxy (e.g. Abadi et al. 2003; Robertson et al. 2004; Gov- 
ernato et al. 2004; 2007). A more serious challenge for the- 
orists is to produce galaxies like M33, as the majority of 
galaxies at these masses do not have significant bulges. For 
a galaxy like M33, with a rotation velocity of 120kms~^ 
and a stellar mass of 1 x 10^" M©, the median n ~ 1.6 and 
45% of blue galaxies have n < 1.5. This means that quasi- 
bulgeless galaxies are the norm in low mass galaxies, which 
means they must form in dark matter haloes with typical 
mass accretion histories. Thus if cosmological simulations 
fail to produce galaxies like M33, this signals a failure of 
the model, and not an unfortunate choice of mass accretion 
history. By contrast, since massive disk dominated galax- 
ies, such as the Milky- Way (which has a bulge fraction of 
~ 0.2 (Widrow, Pym & Dubinski 2008) are rare, they may 
have formed in haloes with a very unusual mass accretion 
histories. This makes simulating the Milky- Way a harder 
problem, because it would take a sample of hundreds of 
Milky- Way mass haloes before one would expect to produce 
a galaxy resembling the Milky- Way. 

4.1 The Wavelength Dependence of the Sersic 
Index 

In Fig. [7] we show the median Sersic index - stellar mass 
relation for the Sersic index measured in each of the 5 SDSS 
pass-bands ugriz. The general trend is for lower Sersic in- 



dices in bluer pass bands. This is true for both red and blue 
galaxies, but the effects are strongest for the blue galaxies. 
This could be attributed to the greater sensitivity of bluer 
pass-bands to extinction and/or variations in the star for- 
mation history as a function of galacto centric radius. The 
lower panels of Fig.[7]show the differences between the Sersic 
index as measured in the u,g,r, and z-bands to that mea- 
sured in the i-band. The largest differences occur around the 
transition mass of 2.5 x lO^^M©. 

These results have implications for studies of galaxy 
structure at high redshift. Currently large surveys with HST 
imaging, such as the Galaxy Evolution and Morphology 
Survey (GEMS) (Rix et al. 2004), the All- wavelength Ex- 
tended Groth strop International Survey (AEGIS) (Davis 
et al. 2007), and the cosmological evolution survey (COS- 
MOS) (Scoville et al. 2007; Koekemoer et al. 2007) are lim- 
ited to the F814W or the F850LP filters (hereafter referred 
to as the J814 and z850 bands). At redshift z = 1 the 1814 
filter shifts into the f/-band, while the z850 filter corresponds 
to B-band. In an attempt to make a fair comparison be- 
tween galaxies at different redshifts, authors often compare 
structural properties in the rest frame B or V bands at low 
and high redshifts. However, given that these bands do not 
trace stellar mass even at low redshift, and that star forma- 
tion rates at z ~ 1 were an order of magnitude higher than 
they are today (e.g. Noeske et al. 2007), it is not clear that 
such a comparison gives an unbiased view of the evolution of 
galaxy structure in stellar mass. To determine the structural 
properties of galaxies at 2; > 1 free of the biases of young 
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Figure 8. Correlations between color and Sersic index. The upper left panel shows color vs Serisc index for all galaxies in our sample, 
color coded by whether they are in the red sequence or blue cloud. This shows a clear correlation between color and Sersic index. The 
upper right panels show the residuals of the color-mass relation versus the residuals of the Sersic mass relation, both for red and blue 
galaxies. The solid lines show the best fit of A color vs A Sersic index, with dashed lines indicating the error on the slope. The correlation 
coefficient, r, for each fit is also shown. This shows that within the blue cloud or red sequence, the scatter in color and Sersic index, at a 
given stellar mass, are uncorrelated. The lower panels show the color-Sersic relations for 6 mass bins. The range in log Mstar /[h^Q Mq] 
for each bin is given in the top left corner of each panel. The correlation coefficient between color and Sersic index is given in the top 
right corner of each panel. The median color and Sersic index for each mass bin is indicated by a cross. This shows that whether a galaxy 
is red or blue does depend on its Sersic index, especially near the transition mass of Mstar ^ 2.5 X 101°Mq. 



stellar populations and dust will require high resolution NIR 
imaging, preferably in more than one band to a provide rest 
frame red pass band and an optical color to correct for stellar 
population and extinction gradients within galaxies. 



4.2 Is Sersic Index Correlated with Color? 

From Fig. [5] we see that higher mass galaxies are both redder 
and more concentrated than lower mass galaxies. Thus glob- 
ally there is a correlation between Sersic index and color: 
galaxies with lower Sersic index are bluer; galaxies with 



higher Sersic index are redder. This correlation is shown in 
the upper left panel of Fig. [S] However, this does not mean 
that there is a causal relation between Sersic index and color. 
The upper right panels of Fig. |S] show the residuals of the 
color-mass relation vs. the residuals of the Sersic-mass rela- 
tion. For both red and blue galaxies the correlation is con- 
sistent with zero, i.e. within the blue cloud or red-sequence, 
the scatter in color, at a given stellar mass, is uncorrelated 
with the scatter in Sersic index. How do we reconcile this 
fact with the strong global correlation between Sersic index 
and color? 
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Figure 9. Concentration of the I-band light versus stellar mass for three feedback models. The black lines show the 16th, and 84th 
percentiles (dashed), and a fit to the median values (solid). The red line shows the median relation for the observations (from Fig.[5]after 
converting from Sersic index to concentration. The dotted lines show, for reference, the concentration corresponding to a Sersic index of 
n = 1,2, and 4. The energy driven feedback model produces disk galaxies that have Sersic indices that are in good agreement with the 
observations. 



The middle and lower panels of FigH] show the color- 
Sersic relation for 6 mass bins from Mstar = lO**** to 
W^^'^Mq. This shows that at a given stellar mass whether 
a galaxy is blue or red does depend on its Sersic index, es- 
pecially near the transition mass of 2.5 X 10^°Mq. Galaxies 
with high Sersic indices are more likely to be on the re- 
sequence than in the blue cloud. This result is consistent 
with Bell (2008), who showed that having a bulge is a re- 
quirement for a central galaxy being on the red-sequence, 
which is a central prediction of AGN feedback scenarios for 
quenching star formation. 



5 COMPARISON BETWEEN MODELS AND 
OBSERVATIONS 

Having established observationally that Sersic index is a 
strong function of stellar mass, we now investigate whether 
our models can reproduce this trend. We consider three 
models with different feedback prescriptions: 1) no feed- 
back; 2) momentum driven feedback; and 3) energy driven 
feedback. The parameters of these models are given in Ta- 
ble [21 For each model we generate Monte Carlo samples 
with halo masses ranging from IQi" < Mvir < 10^^ h'^ Mq, 
corresponding to virial velocities ranging from 31 < Vvir ^ 
310kms~^. We also include log-normal scatter in the halo 
concentration, c, halo spin parameter, A, and halo angu- 
lar momentum shape parameter, a. The scaling relations of 
these models are discussed in more detail in DB08. 

5.1 Concentration Mass Relation 

The upper panels of Fig. [9] shows the concentration of the 
Tband light as a function of stellar mass for these three 
models. Recall that we define the concentration parameter 
as the ratio between the radii enclosing 80 and 20% of the 
I-band light: C82/ = Rsoi/R20i- 

The model without feedback results in a median concen- 
tration of ~ 9, which corresponds to a Sersic index of ~ 3.5. 
The high concentrations are a result of the high galaxy mass 



Table 2. Model Parameters 



Name 


tEFB 


eMFB 


A 


0"lnA 


AC 


I: No Feedback: 


0.0 


0.0 


0.035 


0.35 


N 


II: Momentum: 


0.0 


1.0 


0.035 


0.35 


N 


III: Energy: 


0.25 


0.0 


0.025 


0.35 


N 



fractions, which cause the baryons to dominate the inner 
circular velocity profile. The model with momentum driven 
wind also results high concentrations at high stellar masses 
Mstar — 10^^ Mq, but also shows a trend of decreasing con- 
centration with decreasing stellar mass, in qualitative agree- 
ment with the observations (solid red line). However, the 
observed relation has a sharper transition, as well as both 
lower concentrations at low masses and higher concentra- 
tions at high masses. The energy driven wind model results 
in a much better agreement with the observed relation. It has 
the correct asymptotic values for the concentrations at both 
low and high masses. In particular this model successfully 
produces galaxies with close to exponential surface bright- 
ness profiles at low stellar masses. We emphasize that none 
of these models have been tuned to reproduce the concen- 
tration mass data. The feedback efficiencies of these models 
were chosen so that the models reproduce the zero points of 
the TF and size-stellar mass relations (see DB08). 

5.2 Scatter in Concentrations 

In all three models there is a significant scatter in the con- 
centration index at a given stellar mass. To determine where 
this scatter comes from we show in Fig.[TO]the concentration- 
mass relation for the energy driven feedback model where 
there is only one source of scatter: c. A, or a. This shows that 
the scatter in the I-band concentration-stellar mass relation 
is dominated by the angular momentum shape parameter, 
with the halo spin parameter contributing some scatter, and 
halo concentration resulting in a negligible amount of scat- 
ter. Note also that the scatter in concentrations depends on 
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Figure 10. Concentration-mass relation for the energy feedback model for three different sources of scatter: c (left); A (middle) and a 
(right). The scatter in the concentration mass relation is dominated by scatter in a, the angular momentum shape parameter. The lines 
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Figure 11. Difference between concentration measured in I-band to that measured in stellar mass (black), K-band light (red) and 
B-band light (blue), vs stellar mass. The energy driven feedback model produces disk galaxies that have color gradients that are in j 
agreement with the observations. 



stellar mass, with less scatter at lower masses. This shows 
that feedback is able to, at least partially, erase the initial 
conditions. 

The scatter in the concentration mass relation in the 
models is significantly larger, by at least a factor of two, 
than in the observations. In Fig. [1] we saw that the SS05 
AMD results in larger scatter in disk concentration that the 
BOl AMD. Thus an important question for theorists is to 
determine in detail the range of AMD's of dark matter and 
gas in cosmological simulations. 



may be surprising is that the K-band concentrations are not 
equal to those of the stellar mass. Furthermore, the differ- 
ence between I-band and stellar mass is about the same as 
that between B-band and I-band. This suggests that an ac- 
curate determination of the concentration of the stellar mass 
in disk galaxies requires the color gradients (and hence stel- 
lar mass-to-light gradients) to be taken into account. 



5.3 Differences Between Concentrations in Stellar 
Mass and Light 

Fig llll shows the differences, with respect to the I-band light, 
of the concentration indices measured in the stellar mass 
(black lines), K-band light (red lines) and B-band light (blue 
lines). As might be expected, in all three models the con- 
centrations are highest in stellar mass, and lowest in B-band 
light. As discussed in ii l3.6l the color gradients are a reflection 
of the inside-out formation of stellar disks. However, what 



In the observations, low mass galaxies, Mstar — 10® Mq, 
have small differences between Sersic indices in different pass 
bands. The differences between Sersic indices in different 
pass bands increase with stellar mass, reaching a maximum 
at around the transition mass of 2.5 X lO^°A/0. Qualita- 
tively similar trends of light concentration with stellar mass 
are found with our energy feedback model. The momentum 
feedback and no feedback models, however, have almost no 
mass dependence to the concentrations, in clear conflict with 
the observations. 
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Figure 12. Examples of galaxies with exponential stellar density profiles (solid blue lines) in our energy driven feedback (eEFB = 0.25) 
model. Virial mass increases left to right from Mvir = 10^'^h~^ Mq to Mvir = 10^^h~^ Mq. All models follow the concentration mass 
relation. The halo angular momentum parameters, a and A, for each model are specific at the top of the surface density panels. These 
values are typically within Itr of the values predicted from cosmological simulations. The models produce exponential stellar disks with 
(lower panels) and without (upper panels) adiabatic contraction. The line types and colors are the same as in Fig[T] 



5.4 Examples of Model Exponential Stellar Disks 

From the light concentration - stellar mass plots of Figs|9] 
& [To] it is apparent that galaxies with concentrations con- 
sistent with pure exponential surface brightness profiles are 
produced by our models with feedback, typically within the 
la variations of the halo angular momentum parameters. In 
Fie: ll2l we show examples of such cases. These models have 
exponential stellar disks over the radial range of ~ to ~ 5 
stellar disk scale lengths. Virial mass increases left to right 
from Afvir = lO^°/i"^M0 to Mvh- = W^^h'^MQ. The halo 
angular momentum parameters (q, A) of these models are 
specified above the surface density plots. For comparison, 
the median and la upper scatter in these values in cosmo- 
logical simulations are (0.90,1.15) for a (Sharma & Stein- 
metz 2005) and (0.035, 0.60) for A (Maccio et al. 2007). The 
upper panels are galaxies taken from Model III, while the 
lower panels are taken from an equivalent model but with 
adiabatic contraction of the halo. The models with adiabatic 
contraction also produce pure exponential stellar disks, al- 
beit with slightly higher values of {a and A). This demon- 



strates that the ability of our outflow models to produce 
pure exponential stellar disks is not dependent on the effect 
of adiabatic contraction. 

5.5 Are Outflows Needed to produce Exponential 
Stellar Disks? 

We have shown that a model with SN driven mass outflows 
can reproduce the observed dependence of light concentra- 
tion with stellar mass of spiral galaxies. In particular this 
model can produce pure exponential stellar disks over a wide 
range in halo masses (Mvir = 10^" - lO^^/i^^ Mq). We now 
ask the question of whether this result requires mass out- 
flows, or whether it could be achieved by a model in which 
the baryon fraction was regulated by another mechanism. In 
another words: Is getting the right amount of cold baryons 
in a galaxy sufficient to produce exponential stellar disks? 

We test this hypothesis by running models with the 
same galaxy mass fraction at redshift zero as our exponen- 
tial stellar disk models from Fig 12. We implement this in 
our model by turning off the cooling time calculation, and 
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Figure 13. Feedback modifies the baryon fraction, total angular momentum and angular momentum distribution of the baryons. The 
models on the far left are the Afvir = 10^^h~^ (upper panels) and Mvir = lO^^h"^ TsAq (middle panels) and and A/vir = 10^^h~^ Wq 
(lower panels) models from Fig |12l (which were chosen to have exponential stellar disks). These models have energy driven feedback. The 
models second from the left have been constructed to have the same rrigai and Ag^i as the models on the far left. These models have 
different baryonic and stellar disk density profiles, with higher concentrations and smaller stellar disks in the models without feedback. 
This shows that feedback is modifying the angular momentum distribution. The models second from the right also have mgai fixed to the 
value obtained with feedback, but has Ag^i fixed to the A of the halo. The models on the far right have no feedback, and thus the rTig^i 
and Agai are determined by the efficiency of cooling. These models have too many baryons, which results in disks that are too small, too 
concentrated, and with unrealistic rotation curves. See text for further details. 
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instead we set the amount of baryons that accrete onto the 
disk, at each time step, to be the desired galaxy mass frac- 
tion. The result is a model with our desired galaxy mass 
fraction as well as a galaxy spin parameter equal to the in- 
put halo spin parameter. We run this procedure for three 
halo masses: Mvir = 10^", lO", & IO^^/i'^Mq. The results 
of these test are shown in Fig. [13] 

The models with feedback (and exponential stellar 
disks) are on the far left. Models with the same galaxy mass 
fractions and galaxy spin parameters (set by hand) are in the 
panels second from the left. These models have smaller disks 
with higher concentrations. This demonstrates that outflows 
have modified the distribution of angular momentum of the 
baryons relative to that of the dark matter. The models sec- 
ond from the right have galaxy mass fractions fixed to those 
of the models with feedback, but with galaxy spin parame- 
ters set equal to that of the dark matter. The models on the 
right have no feedback, and thus the mgai and Agai are de- 
termined by the efficiency of cooling, which is close to 100% 
for these halo masses. These models have too many baryons, 
which results in disks that are too small, too concentrated, 
and with unrealistic rotation curves. 

This sequence of models shows that getting the right 
amount of baryons into a galaxy is not sufficient to produce 
exponential stellar disks. Furthermore, getting the right 
amount of baryons and specific angular momentum is also 
not sufficient to produce exponential stellar disks. For the 
models to produce exponential stellar disks they need to get 
the right amount of baryons, specific angular momentum, 
and distribution of specific angular momentum. SN driven 
mass outflows are a mechanism that reduces galaxy mass 
fractions and at the same time naturally modifles the an- 
gular momentum distribution of the baryons in the desired 
way. It remains to be seen whether other mechanisms for 
reducing the baryon fraction of galaxies, can be found that 
also modify the angular momentum of baryons in the desired 
way. 



6 DISCUSSION 

6.1 Are Galaxy Disks Universally Exponential? 

A common assumption in galaxy formation models, and ob- 
servational bulge-to-disk decompositions is that disks are 
exponential. This is not without some motivation. 1) Galax- 
ies that are observed to be bulge-less (i.e. they have no 
spheroid) , invariably have close to exponential light profiles, 
NGC 300 is a example of a galaxy with an exponential stel- 
lar disk over 10 scale lengths (Bland-Hawthorn et al. 2005). 
2) Even if a galaxy is not bulge-less, the disk is often well 
fitted by an exponential over several scale lengths, so it is 
tempting (especially when plotting surface intensity on a 
log-linear plot) to extrapolate the exponential disk to small 
radii, and to associate all of the excess light to a bulge. 3) 
Exponential profiles are also easy to work with analytically: 
mass profiles, and rotation curves can be computed easily. 
However, there are several reasons to question the assump- 
tion that stellar disks are universally exponential. 

From a theoretical point of view, the formation of pure 
exponential disks is a mystery. It has been known since Mes- 
tel (1963) that a uniform sphere in solid body rotation has a 



specific angular momentum distribution close to that of an 
exponential disk. But as shown by e.g. Dalcanton Spergel & 
Summers (1997), when a gas disk with a solid body AMD is 
in centrifugal equilibrium in an dark matter halo, the center 
of the disk is significantly more concentrated than an expo- 
nential, and there is an outer cut off at around 4 disk scale 
lengths. 

Viscous evolution is often cited as a mechanism for ex- 
plaining exponential disks (Olivier, Blumenthal, & Primack 
1991; Ferguson & Clarke 2001; BeU 2002). However, it re- 
quires initial conditions that are less concentrated than ex- 
ponential to start from, which in ACDM is the problem. 
Viscous process always increase the amount of mass at small 
radii, so even if it is a mechanism to explain why stellar disks 
are exponential over a large number of scale lengths, it does 
not explain the existence of bulge-less exponential disks. 

Angular momentum distributions of baryons and dark 
matter in a ACDM cosmology are typically quite different 
from that required to produce exponential disks (Bullock 
et al. 2001b; Sharma & Steinmetz 2005). In particular, they 
tend to have too much low and high angular momentum ma- 
terial. The excess of high angular momentum material may 
not be a serious problem as at larger radii gas dominates 
over stars. In addition gas at large radii can be ionized mak- 
ing it hard to detect. The excess of low angular momentum 
material is considered the more serious problem. 

We have shown that it is possible (through SN driven 
outflows and density dependent star formation) to make disk 
galaxies with realistic stellar surface brightness profiles with 
ACDM initial conditions. In many cases the disk is exponen- 
tial over several scale lengths, but upturns above an expo- 
nential profile at small radii are ubiquitous. When present in 
observed galaxies, these upturns above an exponential, are 
generally considered to be due to a bulge or a pseudo-bulge. 
However, given that there is no theoretical motivation for 
the existence of pure exponential disks, it seems reasonable 
to suppose that in some cases these upturns are actually part 
of the disk. Such a conclusion has been reached before: quot- 
ing from Kormendy & Kennicutt 2004 "We already accept 
Freeman 1970 Type II profiles as canonical disk behaviour, 
even though we can explain it in only a few cases (e.g. Tal- 
bot, Jensen & Dufour 1979). We accept outer cutoffs (e.g. 
van den Kruit & Searle 1981a, b, 1982). Oval disks are only 
piece-wise exponential. Would it be a surprise if disks also 
knew how to deviate above an exponential to small radii?" 

6.2 On the Origin of the Hubble Sequence 

Observationally it is well established that galaxy concentra- 
tion (e.g. Kauffmann et al. 2003); bulge fraction (e.g. Ben- 
son et al. 2007) and Sersic index (e.g. this paper) increase 
with stellar mass. An alternative way to frame these obser- 
vations is that the Hubble sequence of galaxy types is also, 
on average, a mass sequence. Understanding the origin of 
the Hubble sequence is one of the most fundamental prob- 
lems in astrophysics. Here we identify four mechanisms that 
give plausible explanations for the observed trends. All of 
these mechanisms likely occur, but it remains to be deter- 
mined which, if any, process dominates on any particular 
mass scale. 

1) Major Merger Rate: Based on the differences in 
the shape of the halo mass function to that of the stel- 
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lar mass function, Mailer (2008) showed that major galaxy 
mergers are more common in higher mass galaxies. For ex- 
ample, the largest contribution in mass comes from mergers 
with a mass ratio of 1:10 (e.g. Stewart et al. 2008). In a halo 
with mass less than about IO^'^Mq a 1:10 halo mass merger 
ratio will result in a less than 1:10 galaxy mass merger ra- 
tio, whereas in a halo of mass greater than about 10^^ Mq 
a 1:10 halo mass merger ratio will result in a greater than 
1:10 galaxy mass ratio. Since major mergers destroy stel- 
lar disks producing spheroidal remnants (e.g. Barnes 1992; 
Cox et al. 2006), higher mass galaxies are expected to have 
higher bulge fractions. 

2) Gas Fractions: It is observationally known that cold 
gas fractions are lower in higher mass galaxies (McGaugh & 
de Blok 1997; Kannappan 2004). For disk galaxies this is 
likely caused by the combined effects of the Schmidt law for 
star formation (in which star formation is less efficient at 
lower gas densities) and the effects of outflows (which result 
in lower density disks in lower mass galaxies). Theoretically 
it is known that the bulge fraction of the merger remnant 
increases with decreasing gas fraction of the progenitors (e.g. 
Springel & Hernquist 2005; Robertson et al. 2006; Hopkins 
et al. 2009). Thus we expect higher mass galaxies to have 
higher bulge fractions due to major and intermediate mass 
ratio mergers. 

3) Secular Evolution: Secular (i.e. slow relative to the 
dynamical time of the system) processes such as bars, oval 
distortions and spiral structure are known to redistribute 
mass within disk galaxies (see Kormendy & Kennicutt 2004 
for a review). A simple criterion for the onset of disk insta- 
bility is the relative contribution of the disk to the total rota- 
tion velocity within the optical radius of the galaxy (e.g. Ef- 
stathiou, Lake, & Negroponte 1982). Since the contribution 
of the disk to the total rotation velocity increases with stellar 
mass and surface brightness (e.g. Zavala et al. 2003; Pizagno 
et al. 2005; McGaugh 2005; Button et al. 2007), more mas- 
sive disks are, on average, expected to be more prone to disk 
instabilities. Thus more massive galaxies should have larger 
bulge fractions and light concentrations than less massive 
galaxies. 

4) Feedback: In our model the structure of the disk 
is determined by conservation of specific angular momen- 
tum. The mass dependence to the disk structure arises due 
to two effects: 1) the specific angular momentum of galax- 
ies increases with decreasing halo mass, due to increased 
outflow efficiency in lower mass haloes. 2) Feedback pref- 
erentially removes low angular momentum material, which 
directly results in less concentrated disks. This model pro- 
duces an increase in the concentration of the stellar disk 
with stellar mass, in agreement with observations, but with- 
out secular evolution or mergers. It is also worth pointing 
out that in order to produce galaxies with the observed gas 
fractions and standard star formation laws, feedback is re- 
quired (Button & van den Bosch 2008). 



6.3 Future Directions 

There are a number of avenues for future research suggested 
by this paper. 



6. 3. 1 Observations of Disk Galaxy Structure 

The observational analysis of galaxy structure that we 
present in §4 could be greatly improved upon. The first step 
would be two dimensional Sersic fits, which take account of 
the inclination of the disk. However, single Sersic fits do not 
always provide a good match to the data, for example if 
there are multiple components with different axis ratios, or 
if there is a highly concentrated bulge/point source in the 
center of the galaxy (which will drive the Sersic index ar- 
tificially high). An improvement would be multi-component 
fits, such as bulge plus disk (e.g. Simard et al. 2002), or 
bulge plus disk plus bar fits (e.g. Weinzirl et al. 2008). 

Rather than focusing on bulge to disk ratios, which are 
notoriously hard to determine (and depend implicitly on 
how one defines the bulge and disk), an alternative, and per- 
haps more useful approach would be to measure the average 
(and dispersion of) de projected surface brightness profiles 
as a function of fundamental galaxy parameters, such as 
stellar mass, color, and size. Not only would such a study 
highlight the role of fundamental galaxy parameters in de- 
termining galaxy structure, but would also provide strong 
observational constraints for galaxy formation models to be 
tested against. 

Since disk galaxies often have strong color gradients, the 
concentration of the stellar mass may be significantly larger 
than that in optical light. Thus it is desirable to measure disk 
structural parameters in wavelengths that are more sensitive 
to stellar mass, and less sensitive to extinction and young 
stellar populations (although these are useful to study in 
their own right). This means near IR rest frame wavelengths. 

The calculation of global stellar masses is now routine 
both at low and high redshift. The next step in galaxy struc- 
ture studies is to apply these methods either to radial surface 
brightness profiles, or directly to the 2-B surface brightness 
maps. The feasibility and importance of these procedures 
have been demonstrated in spiral galaxy mass model studies 
(e.g. Kranz et al. 2003; Kassin et al. 2006). What is needed 
now is a study of the stellar mass structural properties of a 
statistical sample of galaxies. 

An often ignored component of galaxies is the neutral 
and atomic gas. While the cold gas may only contribute 
~ 10% of the baryonic mass of Milky way mass galaxies, 
due to the atomic gas being typically more extended than 
the stars it is likely to contribute significantly to the angular 
momentum. Furthermore, the gas to stellar mass ratio in- 
creases in lower mass galaxies, and in dwarf galaxies the gas 
is often the dominant component. Thus to achieve a com- 
plete observational picture of galaxy structure, it is essential, 
especially in lower mass galaxies, to include the cold gas. 

6.3.2 Improvements to Treatment of Angular Momentum 
in Semi Analytic Models 

A necessary simplification made in our model was the as- 
sumption that the spin and AMB of the halo was, for a 
given halo, constant with redshift. The validity of this as- 
sumption needs to be tested with cosmological simulations. 

We have shown that a subset of the AMBs found in 
cosmological simulations can produce close to exponential 
disks, once the effects of star formation and feedback are 
taken into account. An important question is whether the 
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merger histories of these haloes are consistent with the sur- 
vival of a bulge-less disk. In order to address this question 
one needs to extend the disk evolution models discussed in 
this paper to full halo merger histories (e.g. Stringer & Ben- 
son 2007). However, in order to treat the build-up of angular 
momentum self consistently we need to know how the AMD 
of the halo is correlated to its merger history. For exam- 
ple we would like to split the AMD into that due to major 
mergers (which will generally destroy disks), and that due to 
minor mergers and smooth accretion (which will generally 
result in growth to the disk). 

6.3.3 Initial Conditions for Secular Evolution/Merger 
Simulations 

Numerical simulations of mergers and secular evolution are 
generally (but not always e.g. Debattista et al. 2004) set 
up assuming the disk is exponential (e.g. Springel & Hern- 
quist 2005; Cox et al. 2006; Debattista et al. 2006; Foyle, 
Courteau, & Thacker 2008). The surface density profiles of 
our models, which are motivated by cosmological initial con- 
ditions, are in general not exponential. If these are used as 
initial conditions, do they result in any testable differences in 
galaxy structure, compared to exponential disk initial condi- 
tions. Furthermore, do secular processes, or mergers preserve 
or modify the stellar concentration-mass relation produced 
by our simple disk evolution models? 



7 SUMMARY 

We have used a disk galaxy evolution model to investigate 
the origin of the surface density profiles of disk galaxies. Our 
model follows the accretion, cooling and ejection of baryonic 
mass inside growing dark matter haloes. Contrary to most 
models of disk galaxy formation, we do not assume that 
the disk is an exponential. In our models the surface den- 
sity profile of the disk is determined by the specific angular 
momentum distribution (AMD) of the cooled baryons and 
centrifugal equilibrium. For the specific angular momentum 
of the hot gas we use fitting functions from cosmological 
simulations. The galaxy mass fractions rrigai are determined 
by the efficiencies of cooling and feedback, and the stellar 
mass-to-light ratios are determined by the star formation 
history and chemical evolution of our models. We have used 
observations of Sersic indices from the SDSS to constrain 
our feedback models and provide observational constraints 
on the frequency of bulge-less galaxies. We summarize our 
results as follows: 

• We have compared two fitting formulae for the halo 
AMD: from Bullock et al. (2001b), and Sharma & Steinmetz 
(2005). In general terms these two formulae give similar re- 
sults for the disk structure, but there are significant differ- 
ences, especially at the tails of the distributions. It would 
be desirable to determine which if either of these AMD's 
provides a better description of the AMD's in cosmologi- 
cal simulations, and whether an alternative fitting function 
provides a better overall fit. 

• Exponential and quasi-exponential stellar disks can be 
produced by our model through a combination of SN-driven 
galactic outflows (which preferentially remove low angular 



momentum material), intrinsic variation in the angular mo- 
mentum distribution of the halo gas, and the inefficiency of 
star formation at large radii. 

• We use observations from the SDSS NYU-VAGC to 
show that the median Sersic index of both red and blue 
galaxies is a strong function of stellar mass. For blue galax- 
ies, low mass (Mstar — 10^ Mq) galaxies have n ~ 1.3, while 
high mass (Mstar — 10^^ Mq) galaxies have n ~ 4, with a 
transition mass of Afstar — 2.5 x IO^^Mq. 

• Our model with energy driven outfiows correctly repro- 
duces the observed relation between Sersic index and stellar 
mass, whereas our models with momentum driven outflows 
and no outflows over predict the Sersic indices in low mass 
galaxies. Thus the exponential density proflle is not a trivial 
result of cosmological initial conditions. 

• The observed fraction of "bulge-less" exponential galax- 
ies is a strong function of stellar mass. For Milky- Way mass 
galaxies (Kot ^ 220 kms"\ Mstar lO^M©) less than 
0.1% of blue galaxies are bulge- less, whereas for M33 mass 
galaxies (Kot ^ 120 km s"\ Mstar ^ lO^^M© bulge-less and 
quasi-bulgess galaxies are more common with ~ 45% of 
blue galaxies having Sersic index n < 1.5. These results 
suggest that the difficulty of hierarchical formation models 
to produce bulge-less Milky- Way mass galaxies (e.g. Abadi 
et al. 2003; Robertson et al. 2004; Okamoto et al. 2005; Gov- 
ernato et al. 2007) is in fact not a problem for ACDM. How- 
ever, the problem of producing galaxies like M33 remains, 
and will provide a key test for hierarchical galaxy formation 
models, and in particular the role of feedback. 

• Contrary to general assumptions, our models suggest 
that galaxy disks do not have a universal exponential density 
profile, especially at small radii. 



ACKNOWLEDGEMENTS 

I thank, Frank van den Bosch, Eric Bell, Stephane Courteau, 
Avishai Dekel, Sandy Faber, Neil Katz, Andrea Maccio, and 
Julio Navarro for useful discussions. I thank the referee for 
providing useful suggestions for improving the manuscript. 
A.A.D. acknowledges support from the National Science 
Foundation Grants AST-0507483, AST-0808133. 

Funding for the Sloan Digital Sky Survey (SDSS) has 
been provided by the Alfred P. Sloan Foundation, the 
Participating Institutions, the National Aeronautics and 
Space Administration, the National Science Foundation, 
the U.S. Department of Energy, the Japanese Monbuka- 
gakusho, and the Max Planck Society. The SDSS Web site 
is http://www.sdss.org/ 

The SDSS is managed by the Astrophysical Research 
Consortium (ARC) for the Participating Institutions. The 
Participating Institutions are The University of Chicago, 
Fermilab, the Institute for Advanced Study, the Japan 
Participation Group, The Johns Hopkins University, Los 
Alamos National Laboratory, the Max-Planck-Institute for 
Astronomy (MPIA), the Max-Planck-Institute for Astro- 
physics (MPA) , New Mexico State University, University of 
Pittsburgh, Princeton University, the United States Naval 
Observatory, and the University of Washington. 



21 



22 Button 



REFERENCES 

Abadi, M. G., Navarro, J. F., Steinmetz, M., & Eke, V. R. 2003, 

ApJ, 591, 499 
Abazajian, K., et al. 2004, AJ, 128, 502 

Baldry, I. K., Glazebrook, K., Brinkmann, J., Ivezic, Z., Lupton, 

R. H., Nichol, R. C, & Szalay, A. S. 2004, ApJ, 600, 681 
Barazza, F. D., Jogee, S., & Marinova, I. 2008, ApJ, 675, 1194 
Barnes, J. E. 1992, ApJ, 393, 484 
Bell, E. F. 2002, ApJ, 581, 1013 

Bell, E. F., Mcintosh, D. H., Katz, N., & Weinberg, M. D. 2003, 

ApJS, 149, 289 
Bell, E. F. 2008, ApJ, 682, 355 

Benson, A. J., Dzanovic, D., Frenk, C. S., & Sharpies, R. 2007, 

MNRAS, 379, 841 
Bland-Hawthorn, J., Vlajic, M., Freeman, K. C, & Draine, B. T. 

2005, ApJ, 629, 239 
Blanton, M. R., et al. 2005, AJ, 129, 2562 
Blitz, L., & Rosolowsky, E. 2004, ApJL, 612, L29 
Blitz, L., & Rosolowsky, E. 2006, ApJ, 650, 933 
Blumenthal, G. R., Faber, S. M., Flores, R., & Primack, J. R. 

1986, ApJ, 301, 27 
Bruzual, G., & Chariot, S. 2003, MNRAS, 344, 1000 
Bullock, J. S., Kolatt, T. S., Sigad, Y., Somerville, R. S., 

Kravtsov, A. V., Klypin, A. A., Primack, J. R., & Dekel, A. 

2001a, MNRAS, 321, 559 
Bullock, J. S., Dekel, A., Kolatt, T. S., Kravtsov, A. V., Klypin, 

A. A., Porciani, C., & Primack, J. R. 2001b, ApJ, 555, 240 

(BOl) 

Ceverino, D., & Klypin, A. 2007, arXiv:astro-ph/0712.3285 
Chabrier, G. 2003, PASP, 115, 763 

Chen, D. N., Jing, Y. P., & Yoshikaw, K. 2003, ApJ, 597, 35 

Conroy, C, & Wechsler, R. H. 2008. larXiv:0805.3346l 

Cox, T. J., Dutta, S. N., Di Matteo, T., Hernquist, L., Hopkins, 

P. F., Robertson, B., & Springel, V. 2006, ApJ, 650, 791 
Crampin, D. J., & Hoyle, F. 1964, ApJ, 140, 99 
Dalcanton, J. J., Spergel, D. N., & Summers, F. J. 1997, ApJ, 

482, 659 

Davis, M., et al. 2007, ApJL, 660, LI 

Debattista, V. P., CaroUo, C. M., Mayer, L., & Moore, B. 2004, 
ApJL, 604, L93 

Debattista, V. P., Mayer, L., CaroUo, C. M., Moore, B., Wadsley, 

J., & Quinn, T. 2006, ApJ, 645, 209 
de Jong, R. S., & Lacey, C. 2000, ApJ, 545, 781 
Dekel, A., & Silk, J. 1986, ApJ, 303, 39 
de Vaucouleurs, G. 1959, Handbuch der Physik, 53, 311 
Dommguez-Tenreiro, R., Tissera, P. B., & Saiz, A. 1998, ApJL, 

508, L123 

Dunkley, J., et al. 2009, ApJS, 180, 306 

Dutton, A. A., van den Bosch, F. C, Dekel, A., & Courteau, S. 

2007, ApJ, 654, 27 

Dutton, A. A., & van den Bosch, F. C. 2008. ar Xiv:0810.4963l 

(DB08) 

Efstathiou, G., Lake, G., & Negroponte, J. 1982, MNRAS, 199, 
1069 

Eke, v., Efstathiou, G., & Wright, L. 2000, MNRAS, 315, L18 
Elmegreen, B. G. 1993, ApJ, 411, 170 
Fall, S. M., & Efstathiou, G. 1980, MNRAS, 193, 189 
Ferguson, A. M. N., & Clarke, C. J. 2001, MNRAS, 325, 781 
Firmani, C, & Avila-Reese, V. 2000, MNRAS, 315, 457 
Foyle, K., Courteau, S., & Thacker, R. J. 2008, MNRAS, 386, 
1821 

Freeman, K. C. 1970, ApJ, 160, 811 

Gao, Y., & Solomon, P. M. 2004, ApJ, 606, 271 

Governato, F., et al. 2004, ApJ, 607, 688 

Governato, F., Willman, B., Mayer, L., Brooks, A., Stinson, G., 
Valenzuela, O., Wadsley, J., & Quinn, T. 2007, MNRAS, 374, 
1479 



Gunn, J. E. 1982, Astrophysical Cosmology Proceedings, 233 
Hammer, F., Puech, M., Chemin, L., Flores, H., & Lehnert, 

M. D. 2007, ApJ, 662, 322 
Hoekstra, H., Hsieh, B. C, Yee, H. K. C, Lin, H., & Gladders, 

M. D. 2005, ApJ, 635, 73 
Hopkins, P.,F., Cox, T.J., Younger, J.D., & Hernquist, L. 2009, 

ApJ, 691, 1168 

Kamionkowski, M., & Liddle, A. R. 2000, Physical Review Let- 
ters, 84, 4525 
Kannappan, S. J. 2004, ApJL, 611, L89 

Kassin, S. A., de Jong, R. S., & Weiner, B. J. 2006, ApJ, 643, 
804 

Kauffmann, G., et al. 2003, MNRAS, 341, 54 

Kaufmann, T., Mayer, L., Wadsley, J., Stadel, J., & Moore, B. 

2007, MNRAS, 375, 53 
Kautsch, S. J., Grebel, E. K., Barazza, F. D., & Gallagher, J. S., 

HI 2006, A&A, 445, 765 
Koekemoer, A. M., et al. 2007, ApJS, 172, 196 
Kormendy, J., & Kennicutt, R. C, Jr. 2004, ARA&A, 42, 603 
Kormendy, J., & Fisher, D. B. 2005, Revista Mexicana de As- 

tronomia y Astrofisica Conference Series, 23, 101 
Kranz, T., Slyz, A., & Rix, H.-W. 2003, ApJ, 586, 143 
MacArthur, L. A., Courteau, S., & Holtzman, J. A. 2003, ApJ, 

582, 689 

Maccio, A. V., Dutton, A. A., van den Bosch, F. C, Moore, B., 
Potter, D., & Stadel, J. 2007, MNRAS, 378, 55 

Maccio, A. V., Dutton, A. A., & van den Bosch, F. C. 2008, 
MNRAS, 391, 1940 

Mailer, A. H., Dekel, A., 2002, MNRAS, 363, 1155 

Mailer, A. H. 2008, ASPC, 396, 251 

Mandelbaum, R., Seljak, U., Kauffmann, G., Hirata, C. M., & 
Brinkmann, J. 2006, MNRAS, 368, 715 

Mayer, L., Governato, F., & Kaufmann, T. 2008, Advanced Sci- 
ence Letters, 1, 7 

McGaugh, S. S., & de Blok, W. J. G. 1997, ApJ, 481, 689 

McGaugh, S. S. 2005, Physical Review Letters, 95, 171302 

Mestel, L. 1963, MNRAS, 126, 553 

Mo, H. J., Mao, S., & White, S. D. M. 1998, MNRAS, 295, 319 
Murray, N., Quataert, E., & Thompson, T. A. 2005, ApJ, 618, 
569 

Navarro, J. F., & Benz, W. 1991, ApJ, 380, 320 
Navarro, J. F., & White, S. D. M. 1994, MNRAS, 267, 401 
Navarro, J. F., Frenk, C. S., & White, S. D. M. 1997, ApJ, 490, 
493 

Navarro, J. F., & Steinmetz, M. 2000, ApJ, 538, 477 
Noeske, K. G., et al. 2007, ApJL, 660, L43 

Okamoto, T., Eke, V. R., Frenk, C. S., & Jenkins, A. 2005, 
MNRAS, 363, 1299 

Olivier, S. S., Primack, J. R., & Blumenthal, G. R. 1991, MN- 
RAS, 252, 102 

Pizagno, J., et al. 2005, ApJ, 633, 844 

Pohlen, M., & Trujillo, I. 2006, A&A, 454, 759 

Rix, H.-W., et al. 2004, ApJS, 152, 163 

Robertson, B., Yoshida, N., Springel, V., & Hernquist, L. 2004, 
ApJ, 606, 32 

Robertson, B., Bullock, J. S., Cox, T. J., Di Matteo, T., Hern- 
quist, L., Springel, V., & Yoshida, N. 2006, ApJ, 645, 986 

Roskar, R., Debattista, V. P., Stinson, G. S., Quinn, T. R., Kauf- 
mann, T., & Wadsley, J. 2008, ApJL, 675, L65 

Schmidt, M. 1959, ApJ, 129, 243 

Scoville, N., et al. 2007, ApJS, 172, 1 

Sharma, S., & Steinmetz, M. 2005, ApJ, 628, 21 (SS05) 

Simard, L., et al. 2002, ApJS, 142, 1 

Sommer-Larsen, J., Gelato, S., & Vedel, H. 1999, ApJ, 519, 501 
Sommer-Larsen, J., & Dolgov, A. 2001, ApJ, 551, 608 
Springel, V., & Hernquist, L. 2005, ApJL, 622, L9 
Stewart, K. R., Bullock, J. S., Wechsler, R. H., Mailer, A. H., & 
Zentner, A. R. 2008, ApJ, 683, 597 



© 2008 RAS, MNRAS 000,[TH23] 



stringer, M. J., & Benson, A. J. 2007, MNRAS, 382, 641 
Sutherland, R. S., & Dopita, M. A. 1993, ApJS, 88, 253 
Talbot, R. J., Jr., Jensen, E. B., & Dufour, R. J. 1979, ApJ, 229, 
91 

Tully, R. B., & Fisher, J. R. 1977, A&A, 54, 661 

van den Bosch, F. C. 2001, MNRAS, 327, 1334 

van den Bosch, F. C. 2002a, MNRAS, 332, 456 

van den Bosch, F. C, Abel, T., Croft, R. A. C, Hernquist, L., 

& White, S. D. M. 2002, ApJ, 576, 21 
van der Kruit, P. C, & Searle, L. 1981, A&A, 95, 105 
van der Kruit, P. C, & Searle, L. 1981, A&A, 95, 116 
van der Kruit, P. C, & Searle, L. 1982, A&A, 110, 61 
Wechsler, R. H., Bullock, J. S., Primack, J. R., Kravtsov, A. V., 

& Dekel, A. 2002, ApJ, 568, 52 
Weil, M. L., Eke, V. R., & Efstathiou, G. 1998, MNRAS, 300, 

773 

Weinzirl, T., Jogee, S., Khochfar, S., Burkert, A., & Kormendy, 

J. 2008, arXiv:astro-ph/0807.0040 
Widrow, L. M., Pym, B., & Dubinski, J. 2008, ApJ, 679, 1239 
White, S. D. M., & Rees, M. J. 1978, MNRAS, 183, 341 
Wong, T., & Blitz, L. 2002, ApJ, 569, 157 

Wu, J., Evans, N. J., II, Gao, Y., Solomon, P. M., Shirley, Y. L., 

& Vanden Bout, P. A. 2005, ApJL, 635, L173 
Yang, X., Mo, H. J., van den Bosch, F. C., Pasquali, A., Li, C., 

& Barden, M. 2007, ApJ, 671, 153 
Younger, J. D., Cox, T. J., Seth, A. C, & Hernquist, L. 2007, 

ApJ, 670, 269 

Zavala, J., Avila-Reese, V., Hernandez-Toledo, H., & Firmani, 
C. 2003, A&A, 412, 633 



On the Origin of Exponential Galaxy Disks 



© 2008 RAS, MNRAS 000.[Tlf23l 



